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Abstract 

Available data on measured observables allow deriving estimates for 
some of the phenomenological parameters that characterize the time- 
evolution and decay of the neutral kaon system based on the hypothesis 
of complete positivity. The present experimental uncertainties are still 
too large to permit a full test of complete positivity. 
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Completely positive evolution maps are the most natural tools for the description 
of the time evolution of open quantum dynamical systems. [1, 2] These systems can be 
modeled as being small subsystems in interaction with suitable large environments. The 
global evolution of the closed compound system is described by an unitary map, while the 
reduced dynamics of the subsystem usually develops dissipation and irreversibility. 

Assuming a weak coupling between subsystem and environment, the reduced time- 
evolution is free from non-linear feedback or memory effects; moreover, it possesses very 
basic and fundamental properties, like forward in time composition (semigroup property), 
probability conservation and entropy increase. 

On the states of the subsystem, conveniently described by density matrices, the re- 
duced evolution is realized by a one-parameter (= time) family of maps, transforming 
density matrices into density matrices and with the additional important property of 
being completely positive. [3-6] This set of transformations forms a so-called dynamical 
semigroup. 

Decay systems can be viewed as specific examples of open quantum systems and 
therefore dynamical semigroups can be effectively used in their description. Indeed, the 
time evolution of the neutral kaon system have been recently studied within this frame- 
work. [7, 8] In practice, assuming a complete positive map being at the basis of the dynamics 
of the kaons results in the addition to the standard Wiesskopf-Wigner evolution of an extra 
linear piece, characterized by six real phenomenological parameters. In the following we 
shall explain in detail how bounds on some of these parameters can be obtained from 
available experimental data. 

As usual, we shall model the evolution and decay of the K°-K° system by means of 
a two-dimensional Hilbert space. [9] A convenient orthonormal basis in this space is given 
by the CP-eigenstates \K\) and \K 2 ): 
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= + , \K 2 ) = — [\K») - . (1) 

As already mentioned, states of a quantum system evolving in time can be suitably de- 
scribed by a density matrix p. This is a positive hermitian operator, i.e. with positive 
eigenvalues, and constant trace (for unitary evolutions). In the case of the kaon system, p 
is a two-dimensional matrix. With respect to the basis (1) it can be written as: 



Pi P3 
Pa P2 



(2) 



where p4 = p%, and * signifies complex conjugation. 

The evolution in time of this matrix can be described in general by an equation of the 
following form: 

^ = -iH p{t) + ip(t) & + L\p]. (3) 

The first two pieces in the r.h.s. constitute the standard Weisskopf-Wigner contribution, 
while L is a linear map that is fully determined by the condition of complete positivity. 
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The effective hamiltonian H (the Weisskopf-Wigner hamiltonian) includes a nonher- 
mitian part, that characterizes the natural width of the states: 

H = M-'-r, (4) 

with M and r hermitian 2x2 matrices. The entries of these matrices can be expressed 
in terms of the complex parameters eg, e^, appearing in the eigenstates of H, 

\K s ) = NsQ s ^J , \K L ) = N L (*^ , (5) 

and the four real parameters, ms, 75 and m^, 7l characterizing the eigenvalues of H: 

As =m s - -is , \h = ra L - , (6) 

where N$ = (1 + |es| 2 ) -1 ^ 2 and N L = (1 + (e^l 2 ) -1 / 2 are normalization factors. It proves 
convenient to use also the following positive combinations: 

Ar = 75 - 7 L , Am = m L - m s , (7) 

corresponding to the differences between decay widths and masses of the states K$ and 
Kl, as well as of the complex quantities: 

r± = T ± zAm , Ar ± = Ar ± 2iAm (8) 

withr = ( 7s + 7L )/2. 

The explicit form of the piece L[p] can be most simply given by expanding p in terms 
of Pauli matrices Uj and the identity a^: p = p^cr^, p, = 0, 1, 2, 3. In this way, the map 
L[p] can be represented by a symmetric 4x4 matrix [Z^„] , acting on the column vector 
with components (po, Pi, P2 5 P3)- It can be parametrized by the six real constants a, 6, c, 
a, and 7: [7, 8] 



a b c 
b a (3 
VO c /3 7 



[V] = " 2 



(9) 



with a, a and 7 non-negative. These parameters are not all independent; to assure the 
complete positivity of the map L, they have to satisfy the following inequalities: 

a < a + 7 , Ab 2 < 7 2 — (a — a) 2 , 

a<a + 7, 4c 2 < a 2 - (a - 7) 2 , (10) 

7 < a + a , 4/3 2 < a 2 - (a - 7) 2 . 

(Notice that from the above inequalities automatically follows that the map L is positive; 
indeed, complete positivity is a stronger request than simple positivity.) 



In absence of the additional piece L[p], the evolution equation (3) would transform 
pure states (i.e. states of the form into pure states, in spite of the fact that 

probability is not conserved, since Tr [p(£)]] < Tr [p(0)]. However, this is due to the presence 
of a non-hermitian part in the effective hamiltonian H and not to a mixing-enhancing 
mechanism that makes p(t) less ordered in time. Instead, loss of quantum coherence shows 
up when the extra piece L[p] is also present: it produces dissipation and possible transitions 
from pure states to mixed states. 

Equations of the general form (3), but with different L, have also been used to describe 
the time evolution of the K°-K° system. [10-14, 15] They are based on the general idea 
that the quantum fluctuations of the gravitational field at Planck's length produces loss 
of phase-coherence. [16, 17] In those discussions, among other more obvious requirements 
like conservation of probability and entropy increase, the linear map L[p\ is constrained to 
be simply positive and not completely positive. 

In the case of the neutral kaons, complete positivity of the time evolution j t , trans- 
forming the initial density matrix p(0) into p(t), amounts to the positivity of the natural 
extension of ^ t to a linear transformation on the states of a larger system consisting of 
the kaon system and another finite-level system of arbitrary dimension. Simple positivity 
guarantees that the eigenvalues of any density matrix for the neutral kaon system remain 
positive. Complete positivity is a stronger property in the sense that the same holds for 
density matrices of the compound system. 

Although at first sight the requirement of complete positivity of the evolution equation 
(3) instead of the much milder simple positivity seems a mere technical request, it has far- 
reaching consequences. In particular, when the additional finite-level system alluded above 
is taken to be another kaon system (a physical situation commonly encountered in the so- 
called ^factories), the complete positivity of the time-evolution (3) assures the absence 
of unphysical effects, like the appearance of negative probabilities, that could occur for 
simple positive dynamics. [18] 

In the following we shall compare the predictions of the completely positive dynamics 
(3), (9) with available experimental data on the neutral kaon system. In principle, this 
allows a complete determination of the values of the phenomenological parameters appear- 
ing in the matrix (9), and therefore a test on the condition of complete positivity, i.e. of 
the inequalities (10). In practice, the experimental data are not yet accurate enough to 
test the conditions (10); nevertheless, estimates on the three parameters c, j3 and 7 will 
be obtained. 

The first step is to solve the equation (3) for the density matrix p with arbitrary initial 
conditions. It is convenient to use perturbation theory, assuming a, 6, c, a, (3 and 7 to 
be small, of the same order of magnitude of egAl? and e^AlT. To first order in the small 
parameters, the time-dependence of the four components of the density matrix p(t) are 
explicitly given by: [8] 

Pl (t) =e-*> t p 1 + (e-^* - e-™*) p 2 

+ (el ~ ^) (e^ - e~^) p 3 + (e L - (e"^ - e"^) p 4 , (11a) 
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p 2 (t)=e-^V 2 + ^(e-^-e-^) p x 



2C 

at! 



Ps(t) 



e _ ( r_+A- 7) t p3 _ 



- U*c + 



2C* 



P 4 (t) 



-r_t _ 

2Am 

r -* _ 

iB* 



)p3 + 



2C* 



-r_t 



) P4 



-7s* 



)P1 



2C* 
AT. 



2C 

at! 



2Am 



-r+t 



-7s* 



) Pi - + 



2C 



-r_t 



(« 



-r+t 



e"^*) p 4 , (H6) 



e"^*) p 2 , (He) 



e"^*) p 2 , (lid) 



where pi, P2, P3 and P4 = p%, are the initial values at t = 0, and the following convenient 
notations have been used: 



A = a + a , B = a — a + 2ib , C = c + i(3 . 



(12) 



Furthermore, from (3) and (9) one can extract the contributions pl and ps that 
correspond to the Kl and Ks neutral kaons: 



PL = 



CL + at: 



~ Ar 



c 



and 



PS = 



Ar H 

2C 

Ar + 



2C* 



_47Jc(^g) e L+Ar 
1 



e% + 



2C* 



s at" 



-3- - 8 
Ar ° 



c 



Ar\ 



4fte 



(*) 



(13a) 



(136) 



In these expressions we have kept only the leading contributions in each entry of the two 
matrices; this approximation will suffice for the subsequent discussions. 

The solutions (11) and (13) can now be used to compute certain observables for the 
K° — K° system; indeed, any physical property of the neutral kaons can be extracted from 
the density matrix p(t) by taking its trace with suitable hermitian operators. 

Useful observables are associated with the decays of the neutral kaons into pion states, 
or into semileptonic states kIv. The amplitudes for the decay of a K° state into and 
7r°7r° final states are usually parametrized as follows, in terms of the s-wave phase-shifts 
Si and the complex coefficients Ai, Bi, i = 1, 2: [19] 



A(K° 



7T + 7r" 



1 



(A + B )e ld » + - / =(A 2 + B 2 )e 



A(K° -> ttV) = (A + B ) e iS ° -V2(A 2 + B 2 ) e iS * , 



(14a) 
(146) 



where the indices 0, 2 refers to the total isospin. The amplitudes for the K° decays are 
obtained from these with the substitutions: Ai — > A* and Bi — > —B*. The imaginary 
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parts of Ai signals direct CP-violation, while a non zero value for Bi will also break CPT 
invariance. 

To construct the operators that describe these two pion decays in the formalism of 
density matrices, one has to pass to the Ki, K 2 basis of Eq.(l). It is convenient to label 
the corresponding decay amplitudes as follows: 



A(K 1 -> 7r°7r°) = X 00 , 



A(K 2 -> TT + TT") = y + _ A(K 1 -> TT + TT") , 

^(K 2 - ttV ) = y 00 A(k x - ttV) . 1 1 " 1 



The complex parameters X and Y can be easily expressed in terms of A^ Bi and 8i of 
(14). Taking the phase convention for which A is real, up to first order in the CP and 
CPT violating parameters, one explicitly finds: 



X, 



X, 



oo 



= V2 
V2 



A + ilm(B 
A + ilm(B ) 



J5o 



Ke{A 2 ) + ilm(B 2 ) 
ne(A 2 ) + iXm(B 2 ) 



J82 



while 



where 



e = 



_ ne(B ) 
e i(s 2 -5o) <jie(Ar 



y o^^™-2e' 



V2 



A 



Im(A 2 ) 

l — , . ; + 



A 

Ke(B 2 ) He(B Q ) 



TZe{A 2 ) TZe(A 2 ) 



Ac 



(16a) 
(16a) 

(17) 

(18) 



is the parameter that uniquely signals CP and CPT violation in the transition amplitudes. 

Using (15), the operators that describe the n + n~ and 7r°7r° final states and include 
direct CP and CPT violations are readily found: [15] 



0+- = \x, 



Y. 
\Y. 



+- 



+- 



Coo — l-^i 



00 1 



1 



Y 



00 



Y 

\Y 



00 



00 1 



(19) 



With the help of these matrices, one can now compute the time-dependent decay rate of 
a neutral kaon into two pions. In the case of charged pions, for a generic initial state, one 
has: 



R+-(t) = 



Tr 



p(t)0 



+- 



Tr 



p(0)O. 



(20) 



For a state which is initially a pure K°, the corresponding density matrix evolves in time 
according to (3), with the initial condition 



p(0)=p K o = \K°)(K \ = 



1 (\ 1 
1 1 



In this case, the observable (20) takes the general form: 

R+-(t) = e~ lst + Rl_e~ lLt + 2e~ rt |?? + _| cos( Am t - </>+_) . 



(21) 



(22) 
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The two-pion decay rate for the K L state can be obtained using (13a): 



2C* ^ 
e L + + Y+- 



+ 



7 

Ar 



c 



AT. 



Ar_ 

while the interference term is determined by the combination 

2C 



(23) 



AT. 



+ Y+- 



V+- = W+- 



i<t>+- 



(24) 



Similar formulas hold for the decay of a K° into neutral pions: one just needs to substitute 
the index H — with 00. 

Notice that in the standard quantum mechanical case, the Kl two-pion decay rate 

is quadratic in the parameters and To get a consistent expression when also 

the term L[p] in (3) is present, it is necessary to take into account terms up to second order 
in all small quantities, as in (23). This is the only case among the discussed observables 
for which second order terms are needed. 

The expressions (23) and (24) can be directly compared with the experimental data. 
However, some care needs to be used in the choice of the proper measured results. In fact, 
in fitting their data on the two pion decay, the various experimental groups use the standard 

quantum mechanical formula for (t), obtained from (22) by setting -R+_ = \rj^ | 2 . 

However, as discussed in [15], different experimental setups measured the time dependence 

of -R-i (t) focusing either into the interference region or the long-time tail of the decay 

distribution. In particular, due to a flux normalization, the results of the CERN- Heidelberg 
Collaboration [20] are insensitive to the interference effects; therefore their data can be 
regarded as a determination of R+_: R+- = (5.29 ±0.16) x 10 -6 . On the other hand, the 
attention of the CPLEAR Collaboration [21] is concentrated on the interference region, 

so that their results can be taken to determine both the modulus \rj-\ | and the phase 

0_l of ry-i Their most recent preliminary values for these two quantities read [22]: 

|?7 + _| = (2.316 ± 0.039) x 10" 3 , + _ = (43.5 ± 0.8)°. 

Other interesting observables, directly measured in experiments, are the asymmetries 
associated with the decay into the final state / of an initial K° as compared to the cor- 
responding decay into / of an initial K°. Using the same notations as before, all these 
asymmetries have the general form 



A(t) 



Tr 


PKo(t)O f 




PK°(t)O f 


Tr 


PKo(t)Of 


+ 


PK°{t)Of 



(25) 



where p^o(t) and pK°(t) are the solutions of (3) with the initial conditions of having a pure 
K° and pure K° at t = 0, respectively. These correspond to the initial density matrices: 



Pko = \K°)(K°\ 



1 

-1 



(26) 



and (21) above. 
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One of the simplest asymmetries that can be computed is the one that involves tv + tv~tv° 

in the final state. The operator o responsible for this decay can be obtained following 

steps similar to the ones that led to the expressions of and Coo above. Parameterizing 

the decay amplitudes as: 



A(K 2 



-7T ) = X + _ , A(K ± 



one easily finds: [15] 



O+-0 = \X 



+-o| 



7T + 7r' 



-7T ) = y + _ A(K 2 



7T + 7r" 



+ -0| 



Y 



+-o 



+- 



-7T ), (27) 



(28) 



-o 1 

Inserting this in (25) and using the solutions (11), one obtains, to first order in the small 
quantities, the following expression for the asymmetry integrated over the entire Dalitz 
plot: 

2C 



A+_ (t) =2TZe(e s - 



Ar_ 



2e -Art/2 7^ e ( ??+ _ ) cos(Amt) - Jm(?7 + _o) sin(Amt)] , 



where 



V+-o = es - 



2C 
AfT 



+ Y+- 



(29) 



(30) 



It is precisely the real and imaginary part of this parameter, appearing in the interference 

term of o(t), that is usually fitted with the experimental data. The best available 

determination of rj^ o comes from the CPLEAR Collaboration [22]. In particular, one 

finds: 7£e(?7_| ) = ( — 4 ± 8) x 10~ 3 . (The value of lm(r] + -o) will not be needed for the 

considerations that follow.) 

Other interesting asymmetries involve decays into semileptonic states. The amplitudes 
for the decay of a K° or a K° state into ir~£ + u and tt + £~u are usually parametrized by 
three complex constants x, y and z as follows: [23] 



A(K° 

a(kv 

A(K° 
A(K® 



e+u) =M(l-y) 
ti-£+v) =xA(K° -> 



TV 



TV 



tv~£ + u) , 



(31a) 
(316) 
(31c) 
(31d) 



where M is a common factor. (In Ref.[23] x = z* is used instead of z.) The AS = AQ 
rule would forbid the decays K° — > tv + £~u and K° — > tv~£ + u, so that the parameters x 
and z measure the violations of this rule. Instead, CPT-invariance would require y = 0. 

From the parametrization in (31), one can deduce the decay amplitudes from the 
states (1) of definite CP, and therefore the following expressions for the two operators 
describing the semileptonic decays: 



2 1 
\M\ 2 



y\ 



1 + 2/1 



|l + x| 2 (l + x*)(l-x) 
(l + x)(l-x*) \l-x\ 2 

\z + l\ 2 (z* + l)(z-l) 
(z+l)(z*-l) \z-l\ 2 



(32a) 
(326) 
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The parameters x, y, and z are expected to be very small, and the available experimental 
determinations support this theoretical hypothesis. Therefore, in the following we shall 
treat these constants as the other small parameters in the theory and compute the various 
semileptonic observables keeping only the first order terms in all these parameters. 

With the help of (32), one can compute the following two semileptonic asymmetries, 
for which experimental data are available: 



At(t) = 



Tr 


PKo(t)(O e --O e+ ) 


-Tr 


p K o(t)(o e --O e+ } 


Tr 




+ Tr 


p K0 {t) (o e - 



Tr 


PKo(t)(O e - +0 £+ ) 


-Tr 


pK°(t)(o e - 


Tr 


PKo(t)(Ot-+Oi+) 


+ Tr 


p K o(t)(o e -+O e+ ^ 



(33a) 



(336) 



The first observable is sensible to the K° <-> K° oscillation frequency, while the second one 
(also called A 2 (t)) is the decay rate asymmetry between K° and K° in any semileptonic 
final state. Using the solutions (11) of the kaon evolution equation, one finds the following 
explicit expressions, valid to first order in all small parameters: 



A A m(t) 



N Am (t) 
D(t) 



Mt) 



N e (t) 
D{t) 



(34) 



where the numerators are given by 



N A m(t) = 2e" rt | [1 - (A-i)t] cos(Amt) 



+ 



Tie 



B 



N e (t) = 2Ke ^e L - 
-2e~ rt 
i Ini ( ( /. 



2C* 
AT. 

2C* 



Am 

e~ lst + 21Ze [ e s - 
2C* 



Xm(x + z) 
2C 



Ar! +ne { es Ar_ 



Ar_ 

2C 



sin(Amt) 

-lLt 

cos(Amt) 



(35a) 



Ar. 



Tm e s 



2C_ 
Afl 



Xm(x — z) 



sin(Amt) \ , (356) 



while the common denominator reads 



D(t) 



l-^+ne(x + z) 



e -jst + 



l + ^-ne(x + z) 



(35c) 



(Notice that the expression in (35a) is valid only for "small" times, i.e. as long as 
e -( A -T)* ~ 1 - (A - 

The asymmetry A Am (t) can be used to determine the parameters a and a, via the 
combinations A = a + a and TZe(B) = a — a, once the AS = AQ rule is assumed and 7 has 
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been fixed using other observables. On the contrary, notice that, within our approximation, 
the long-time limit of Ag(t)/2 gives the constant 2TZe(es — 2C/Ar_), which is independent 
from x, y and z. Therefore, this asymmetry allows a determination of this combination, 

independently from the AS = AQ rule, and in alternative to o, see (29). 

Another useful observable involving semileptonic decays is the so called CP violating 
charge asymmetry: 



S(t) = 



Tr 



p(t) (O e+ -O e -) 



Tr 



p(t) (O e+ +O e -) 



(36) 



For a state which is initially a pure K°, it can be explicitly written as: 



5 S e~ lst + 5 L e~ lLt + 2 e~ Tt cos(Am t) 

\ > ~ e -lst + e -7Lt 



(37) 



where 



2C 

5 s = 2ne[e s + ^r- ) - Ke(2y + z - x) , (38a) 
2C* 

d L =2Tle[t L + -^) -Ke(2y + x-z) , (386) 



are the charge asymmetries for the Ks and Kl states, respectively. For sake of simplicity, 
in the expression (37) we have kept only the dominant pieces in each term. Furthermore, 
notice that the two expressions (38) can also be obtained by using directly the expressions 
(13) for ps and pl- 

Although not immediately connected with violations of a specific discrete symmetry, 
by comparing the time-dependent asymmetry 5(t) with the experimental data one can 
in principle measure at the same time the charge asymmetries for both the Ks and Kl- 
Unfortunately, the above expression for 5(t) has not yet been fitted with the experimental 
data, so that one needs to look for separate determinations of 5 s and Sl- This last quantity 
is very well known from the fixed target experiments, and for its value one can take the 
world average given in [24]: 5 L = (3.27 ± 0.12) x 10" 3 

On the contrary, no direct measurement of 5s has been performed yet: the available 
estimates on this quantity are deduced from other observables, e.g. the two semileptonic 
asymmetries At and Aqpt obtained from the general formula (25) when the final state / 
is n + £~u and n~£ + u, respectively. The CPLEAR Collaboration has recently reported a 
value for the combination (At + Acpt)/%, which in ordinary quantum mechanics would 
give precisely an estimate for 5s [22]. Although this is no longer strictly true when the 
generalized evolution equation (3) is used to describe the decay of the neutral kaon system 
(see [8]), in lacking of more complete and precise results we shall take the CPLEAR 
value for (At + Aqpt) /2 as an experimental determination of the Ks charge asymmetry: 
5 S = (3.2 ±1.8) x lO- 3 . 

We are now ready to derive estimates on some of the parameters appearing in the 
dissipative piece L[p] of the evolution equation (3). In this analysis, we shall adopt some 
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simplifying working assumptions. First, we shall ignore the effects of e' in the two-pion 
decay rate (22), since one expects: e'/e^ ~ 10 -4 . Thus, the formulas presented for the 
charged pion decay apply equally well to the neutron pion decay. Also, we shall concentrate 
on possible CPT violations in the time evolution of the neutral kaon system as given by 
the dissipative part (9) and therefore ignore the CPT non preserving contributions in the 

various decay amplitudes; this amounts to set the parameters y and Bq, hence , to 

zero. Furthermore, we shall assume the validity of the AS = AQ rule in the semileptonic 
decays; this choice is supported by theoretical arguments that predict the parameters x 
and z to be neglegibly small. [19] 

From the two and three-pion observables (24), (30) and the two semileptonic charge 
asymmetries (38), one has: 



Tie \e s + 
He (e s - 



2C 

at! 

2C 

atZ 



2 ' 
= ^e(?7 + _o) 



Tie + 
He ( e L - 



2C 

2C 
Af7 



- 8 Jl 
~ 2 ' 

= \v+- \ cos((^) + _) 



(39) 



These relations allow extracting He(es), IZe(eL) and 

1 + r 2 



c 
AT 

A 

Ar 



16 
l + r 2 
16t 



5 s + S L -2 r] + _ cos((p + -) - 2 Tle(r] + _ ) 
6 L - 5s -2 r] + - cos((p + -) + 2 TZe(i] + - ) 



(40a) 
(406) 



where r = 2Am/Ar is the tangent of the so called superweak angle. From the previously 
quoted experimental values of the quantities on the r.h.s. of (38), and taking into account 
that r is very well known, [24] r = 0.9502 ± 0.0034, one obtains: 



c = (1.0 ± 1.4) x 10" 17 GeV , 



(3= (-1.0 ±1.5) x 10" 17 GeV 



(41) 



From these results, one can now estimate the value of the parameter 7, using (23) and 
the determination of the Kl two-pion decay rate from the CERN-Heidelberg Collaboration: 



7 _ R L 

Ar +" 
+ 



(l + r 2 )Ar 



1 + r 2 Ar 2 

[He(e L ) ((r 2 - l)c - 2t0) - lm{e L ) ((r 2 - l)(3 + 2rc)] ; 



(42) 



Xm(ei) can be obtained by taking the imaginary part of (24): Xm(ex,) 
I77+-I sin(0 + _) + 2(rc - /?)/[(! + r 2 )Ar]. In this way, one gets: 



7 = (0.1 ±30.2) x 10" 20 GeV . 



(43) 



Unfortunately, the errors propagation produces in this case a large uncertainty. 
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Estimates for the parameters a and a can be obtained fitting the asymmetry (33a) 
with the available experimental data. Assuming again x = y = z = 0, taking 7 as in 
(43) and for Am the average of [24], we have performed a preliminary % 2 fit of A/± m 
with the experimental data of the CPLEAR Collaboration published in [25]. The result, 
rather independent from the actual value of 7, gives: a = (2.3 ± 2.9) x 10~ 17 GeV and 
a = (2.1 ± 2.9) x 10~ 17 GeV. These values should be taken as indicative; a more accurate 
fit along the lines discussed in [14] should provide more precise determinations. 

The conclusion that can be drawn from this analysis is that the precision of the present 
experimental data on the neutral kaon system is not high enough to allow a meaningful 
test of the inequalities (10) and therefore of the hypothesis of complete positivity. Nev- 
ertheless, as soon as more complete and precise data become available, it will be possible 
to distinguish between simply positive and completely positive time-evolution maps in an 
actual experimental setup. 

In this respect, particularly promising are the planned experiments on correlated neu- 
tral kaons at </>- factories. In fact, entangled kaon states might develop negative probabilities 
under simply positive time evolutions. Since this is forbidden by complete positivity, very 
different behaviours for entangled kaons are predicted in the two cases, with effects that 
might be experimentally detectable [26] . In view of the wide use of dynamical semigroups 
in the description of open systems [1, 2], these checks of complete positivity might have 
consequences that go beyond the specific situation of the neutral kaons. 
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